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We study persistence in one-dimensional ferromagnetic and 
anti-ferromagnetic nearest-neighbor Ising models with paral- 
lel dynamics. The probability P(t) that a given spin has not 
flipped up to time t, when the system evolves from an initial 
random configuration, decays as P (t) ~ l/t Bp with 8 P ~ 0.75 
numerically. A mapping to the dynamics of two decoupled 
A + A — > models yields 8 P = 3/4 exactly. A finite size scal- 
ing analysis clarifies the nature of dynamical scaling in the 
distribution of persistent sites obtained under this dynamics. 
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Recent work on the phenomenon of "persistence" has 
revealed a novel (and in many respects unexpected) way 
by which an interacting many-body system, evolving in 
time, retains memory of its initial state O. Consider a 
one-dimensional ferromagnetic Ising model with nearest 
neighbor interactions, quenched from an initially random 
(infinite temperature) configuration and allowed to relax 
to its global minimum energy (zero temperature) config- 
uration, a state with either all spins up or all spins down. 
Suppose the dynamics is serial, with an attempt to up- 
date a single spin being performed at each time step. Fix 
one spin and ask: What is the probability that this spin 
has not flipped up to time tl This quantity, the persis- 
tence probability, was first found numerically to decay 
as 



P® ~ i 



(1) 



with 9 S a new, non-trivial exponent. (The subscript 
refers to serial dynamics.) The numerical results sug- 
gested 6 S ~ 0.37 0; a later analytic tour de force de- 
rived 9 S = 3/8 exactly |3[ . Later studies have established 
rigorously that the persistence exponent is a new expo- 
nent characterizing the dynamics; it cannot be related 
to either the static exponents v and r\ or the dynamical 
exponent z|ij. 

This paper discusses the problem of persistence in the 
context of one-dimensional Ising models with nearest 
neighbour interactions evolving under parallel dynamics. 
We consider both ferromagnetic and anti-ferromagnetic 
interactions, where the Hamiltonian has the form 



H = —J} (JjQ-i+l. 



(2) 



Without loss of generality, we take J = 1 for ferromag- 
netic interactions and J = — 1 for the anti- ferromagnetic 
case. Each spin can take values +1 (up) or -1 (down). 
The zero temperature dynamics evolves a configuration 
{cr{t)} at time t to a configuration {cr(t+ 1)} at time t+1 
through the following simple rule: For ferromagnetic in- 
teractions, each spin at time t + 1 assumes the value of 
one of its neighboring spins at time t, chosen from right 
or left with equal probability. For anti-ferromagnetic in- 
teractions, the above rule is modified in the following 
way: At time t + 1, assign the value of each spin to the 
negative of the value of one of its neighboring spins, cho- 
sen from right or left with equal probability. Each such 
step in time constitutes a single Monte Carlo step. The 
parallel nature of the dynamics follows from the fact that 
all spins are updated together. 

We have simulated parallel dynamics using the above 
rules on Ising systems of linear size L = 10 2 — 10 6 sites 
and for times t < 10 5 , applying periodic boundary con- 
ditions. We average over a fairly large number of initial 
conditions, typically 10 2 — 10 3 for the smaller lattices 
(L < 10 4 ), starting from configurations in which each 
spin is independently assigned a value 1 or — 1 with equal 
probability. We compute the standard persistence prob- 
ability P(L,t), defined as the probability that the spin 
at a given site in a system of size L has not flipped up 
to time t, averaged over all sites and over an ensemble of 
initial conditions. For L — > oo (in practice for t <C L z ), 
P(L,t) -> P(t). 

Fig. 1 shows the persistence probability P(t) for a fer- 
romagnetic Ising system of linear size L = 10 6 , evolving 
under parallel dynamics. P(t) exhibits a power law tail 
with an exponent p ~ 0.75. The behavior is identical for 
anti-ferromagnetic interactions. For comparison, the cor- 
responding plot for serial dynamics is shown on the same 
figure; the exponent, as advertised, is 6 S = 3/8 to within 
numerical resolution. The exponent S — 3/8 is also ob- 
tained for anti-ferromagnetic interactions under serial dy- 
namics. It is thus natural to guess that 9 P = 3/4 = 26 s , 
and that 9 P as well as 9 S remain unaltered when ferro- 
magnetic interactions are replaced by anti-ferromagnetic 



ones. 




FIG. 1. The persistence probability P(t) in a Id Ising 
model plotted against time t in a logarithmic scale. The lower 
curve (triangles) is for parallel dynamics while the upper one 
(circles) is for serial dynamics. The solid and the dotted lines 
fitted to these curves have slopes 0.75 and 0.375 respectively. 
The slopes remain the same if ferromagnetic interactions are 
replaced by antiferromagnetic interactions. 

Coarsening in ferromagnetic Ising models with serial 
dynamics occurs through the motion and annihilation of 
domain walls separating spins with different orientations. 
Spins located at domain walls have no preference for ei- 
ther orientation. A given spin flips if it is crossed by a 
domain wall. Thus, persistence in this context is equiva- 
lent to the probability that a specified site has not been 
crossed by a domain wall up to time t. If the domain walls 
in this problem are interpreted as particles of type A, the 
following simple reaction-diffusion scheme describes the 
motion and annihilation of domain walls: A + A — » 0, 
with particles diffusing at each time step and annihilat- 
ing on contact 0]. 

The problem of persistence with serial dynamics then 
translates simply into the following: Given a chosen site 
at time t — and an initial configuration of A particles, 
corresponding to domain walls in the initial configura- 
tion, what is the probability that an A particle has not 
crossed that site up to time i? Such a redefinition of 
the problem recasts the question in terms of the essential 
ingredients of the dynamics, the motion and interaction 
of domain walls. It is natural to look for the analog of 
such domain walls in the case of parallel dynamics to 
understand the conjectured relationship 6 P = 3/4 = 29 s . 

Consider first the ferromagnetic case and divide config- 
urations of spins into the following categories: unstable 



spins, implying that they will definitely flip in the next 
time step, stable spins implying that they will not flip 
in the next time step, and "zero-field" spins, which may 
or may not flip, with either possibility occurring with 
probability 1/2. An unstable spin <Ji at site i has both 
neighbors in the state — (Tj, while stable spins point in 
the same direction as their neighbors. Zero-field sites, on 
the other hand, have one neighbor pointing up while the 
other points down. As a consequence, flipping the spin 
at a zero-field site costs no energy. 

These zero-field sites are the analogs, for parallel dy- 
namics, of domain walls in the ferromagnetic Ising model 
with serial dynamics, in a sense which we make precise 
below. Stable sites belong to domains which are either all 
up or all down. A spin in the bulk of such a domain can 
only increase its energy if it flips; it is thus not updated 
through the zero-temperature dynamics. Unstable sites, 
for the case of ferromagnetic interactions, are associated 
with an anti-ferromagnetic arrangement of spins of the 
form: . . . 1010101010 . . . , where the notation "1" indi- 
cates an up spin and "0" indicates a down spin. Note 
that all sites interior to such a region will flip in the 
next time step. Within a region of unstable sites, spin 
histories follow a two-cycle; each spin flips once in each 
time-step. Sites within such regions cannot contribute to 
persistence, for they cannot be persistent beyond a single 
time step. It is obvious that persistence of spin config- 
urations at late times can only be associated with spins 
which lie deep within stable regions. 

It is useful for the ensuing discussion to divide the one 
dimensional lattice into two inter-penetrating sublattices 
A and B. For example, we may take all even-numbered 
sites as forming the A sublattice and odd-numbered sites 
as constituting the B sublattice. The state of sublattice 
A(B) at any time t is determined by the state of sublattice 
B(A) at t — 1. The initial states of sublattices A and B 
are uncorrelated with each other. 

We will argue below that the persistence of a site i on 
a given sublattice at times t = 1, 3, 5 . . . (2m + 1) is de- 
termined by the configuration of zero-field sites on the 
other sublattice, while at times t = 2,4, 6... 2m, it is 
determined by the configuration of zero-field sites on its 
own sublattice. With this in mind, it is useful to distin- 
guish zero-field sites on sublattice A from those on sub- 
lattice B. Call every zero-field site on the A sublattice an 
A particle and every zero-field site on the B sublattice, 
a B particle, at time t = 0. Fig. 2 illustrates a sequence 
of configurations at succeeding instants of time through 
which an initial state evolves. The letters A and B in- 
dicate the positions of zero-field sites on the A and B 
sublattices. 

Inspection of Fig. 2 leads to the following understand- 
ing: Zero-field sites are located at (i) the interface be- 
tween two stable regions, in which case this interface is 
composed of two adjacent zero-field sites - this is the case 
for the model with serial dynamics in which case A parti- 
cles represent the bound pair of zero-field sites, (ii) the in- 
terface between a stable and an unstable region, in which 



case there is a single zero-field site and (iii) the interface 
between two unstable regions, in which case there is again 
a pair of adjacent zero-field sites. The motion of zero-field 
sites corresponds to the expansion/contraction of one or 
the other domain it separates. Zero-field sites can also 
annihilate, leading to the coalescence and shrinkage of 
domains, as the system coarsens under the dynamics. 



FIG. 2. Time evolution of the zero-field sites on A and B 
sublattices: the configuration at the earliest time corresponds 
to the bottom row and each successive row is a later time do 
not react with each other. Only particles of the same kind, 
i.e. A or B, can annihilate. 

It is thus obvious that the non-trivial dynamics asso- 
ciated with the persistence phenomenon can only be as- 
sociated with the zero-field sites, for such sites constitute 
the boundaries between unstable and stable regions. It 
is also obvious that sites which are persistent till time t 
can only be associated with regions which are stable upto 
time t, i.e. regions which are stable at t = but which are 
not crossed by a zero-field site upto time t. The problem 
of persistence in the case of parallel dynamics can thus 
be discussed precisely as in the serial case, with the dy- 
namics of the simple domain walls in the serial case being 
replaced by a somewhat more complex dynamics in the 
parallel case, which we deal with below. 

We now establish the following crucial ingredient of 
this model which enables the exact calculation of 9 P , us- 
ing the known result for 9 S and the mapping onto the 
equivalent reaction diffusion problem: Only particles of 
the same type (A or B) can annihilate each other. They 
are transparent to particles of the other type. Thus, the 
problem of the motion of zero-field sites in this model 
reduces to the problem of two decoupled A + A — > and 
B + B — > models, with the initial number of A's and 



B's being set by the initial conditions. This is clearly ev- 
ident from Fig. 2 if we look at the configurations at even 
or odd time steps and at even or odd lattice sites. Our 
use of periodic boundary conditions implies that N(A) 
and N(B) are always even, where N(A), N(B) represent 
the number of particles of type A or B at any time. 

The single conservation law of the A + A — > 
model, the conservation of particle number N(A) mod 
2, is replaced by two conservation laws in this model: 
mod(N(A),2) = and mod(N(B),2) = 0, where N(A) 
and N(B) are the numbers of A and B particles in any 
configuration of the model. It is useful to note the fol- 
lowing: while the properties conventionally computed for 
the reaction diffusion scheme A + A — > refers to sequen- 
tial dynamics, these results generalize trivially to parallel 
dynamics. This is a consequence of the fact that particle 
moves are independent of each other in this model. Only 
the reaction itself, in which two A particles on the same 
site annihilate each other, has different interpretations in 
parallel and serial dynamics, but this ambiguity is easily 
removed through a simple redefinition of the time and 
length-scales. 

For the A + A — » model, the density of A particles 
decreases as B 



N(A) ~ l/t 1 ' 2 



(3) 



In the context of the Ising model with parallel dynam- 
ics, this result implies that N(A) and N(B) as defined 
above and at times t, t + 2, t + 4 . . . are conserved modulo 
2 and decay separately as l/t 1 ' 2 , with a prefactor which 
depends on the initial concentration. This result accords 
with results obtained for this model by PrivmanB; we 
have checked this numerically as well. Also note the fol- 
lowing: If we consider only the persistence of sites on 
a given sublattice, at only even or odd time steps i.e. 

t = 0,2, ,4. . . .2m or t = 1,3, 5 2m + l, this will decay 

as P(t) ~ 1/t 3 / 8 , reflecting the fact that this dynamics 
maps exactly onto the A + A — > dynamics. This is 
consistent with our numerics. 

We now use these results to argue the following. Con- 
sider, for concreteness, the persistence of a spin on the 
A sublattice at some time t after the quench from in- 
finite temperature. The persistence probability is then 
simply the probability that the site in question has not 
been crossed by an A particle upto time t or a B particle 
upto time t — 1. Since these probabilities are indepen- 
dent, (the dynamics of A and B particles decouples), this 
joint probability is simply the product of the independent 
probabilities that the site is persistent with respect to the 
motion of both A and B particles, implying 



P(t) 



1 



£3/8 ( t _ !)3/8 



t 3/4 



(4) 



yielding the persistent exponent for parallel dynamics 
Opar = 3/4 exactly, consistent with the numerical data. 

To test the validity of the mapping onto the two 
non-interacting species of particles (A and B) out- 



lined above, we have simulated the associated reaction- 
diffusion model independently and computed, numeri- 
cally, the analog of the persistence probability for the 
Ising case. This is done by computing the probability 
that a given site is crossed by a particle of neither type 
upto time t; the exponent obtained numerically tallies 
precisely with our result above. 

How do these results generalize to the anti- 
ferromagnetic Ising model with parallel dynamics? These 
results are unaltered as a consequence of the following 
simple mapping of configurations: Replacing all spins on 
one sublattice, say the A sublattice, through {a a} — > 
—{a a}, changes the sign of the exchange interaction J, 
mapping the problem with the new variables into the fer- 
romagnetic problem. This gauge symmetry relates con- 
figurations pairwise; every update for the ferromagnetic 
case is an allowed update for the anti-ferromagnet with 
the same weight. Thus none of the conclusions here are 
altered and the persistence exponent is independent of 
the sign of the exchange interaction J. 

Recently, there has been considerable interest in the 
spatial scaling properties of persistence 0. Numerical 
work on the 1-dimensional A + A — > model, which de- 
scribes the Ising model with serial dynamics, shows the 
existence of a non-trivial fractal structure in the spatial 
distribution of persistent sites at long times. These re- 
sults can be recast in terms of a dynamical scaling form 
for P(L,t)§, 



version of the model, a result corroborated by the work 
of Bray and O'Donoghue || and Manoj and Ray J7J. 



P{L,t) = lT^f(t/L'), 



(5) 



where z is the dynamic exponent and f(x) ~ x~° s for 
x << 1 while it is constant for large x. One consequence 
of this form is that persistent sites at long times and for 
length scales £ << t 1 / 2 constitute a fractal with fractal 
dimension df = d — z9 s — 0.25. We have used z = 2, 
valid for A + A — > dynamics. 

Does such structure exist for the parallel version of 
Ising persistence? Our data for P{L,t) are consistent 
with Eq.rtq), with 9 S replaced by 9 P and z ~ 2, as shown 
in the scaling plot of Fig. 3. This illustrates the valid- 
ity of the dynamical scaling ansatz for persistence under 
parallel dynamics. (Data collapse here is, however, infe- 
rior in comparison to the serial case). Using the result 
of the previous paragraph, this would indicate a fractal 
dimension of —0.5, a-priori an unphysical result. This 
result can be attributed to the fact that we are looking 
at the persistence of a site under two independent pro- 
cesses. Consider the set of sites persistent with respect 
to the motion of A and B particles separately; each will 
form a fractal with the same fractal dimension df. The 
intersection of these two fractals represents those sites 
persistent with respect to the motion of both A and B 
particles. Assuming that A and B particles are initially 
uncorrelated, the dimension of the intersection set is then 
2df — d = —0.5, as above. We conclude that persistent 
sites in the parallel dynamics version of the Ising model 
do not exhibit spatial scaling of the type seen in the serial 
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FIG. 3. Plot of P(L,t)L z0 " vs. t/L z , where P(L,t) is the 
configuration averaged density of persistent sites in a system 
of size L at time t, Z — 2 and 9 P — 3/4, illustrating the 
validity of the dynamical scaling ansatz. The inset shows the 
total number of persistent sites N P (L) = LP(L,t = oo) left 
in the system as t — > oo, plotted against the system size L on 
a logarithmic scale, both for parallel (circles) and sequential 
(squares) dynamics. The straight lines fitted to these points 
have slopes -0.5 and 0.25 respectively (see text for discussion). 

One implication of this result is that both the average 
number and average density of persistent sites in a system 
of size L should decay with L for parallel dynamics, in 
contrast to the serial case. We have verified this numeri- 
cally; see the inset to Fig. 3. (In contrast, the mass of a 
truly fractal object increases with scale while its density 
decreases.) Thus, a large system has no persistent sites, 
at sufficiently long times, for most initial conditions. The 
power-law tail of P(L, t) is then associated with an ever- 
smaller fraction, as L increases, of initial states, whose 
associated persistent sites survive for longer and longer 
times. 

In conclusion, we have studied persistence in one- 
dimensional Ising models with parallel dynamics. We 
have obtained the persistence exponent 9 P = 3/4 exactly 
and studied the spatio-temporal correlations of persis- 
tent sites. These results relied on an exact mapping to 
the dynamics of two decoupled A + A — > models. It 
would be interesting to see if similar arguments exist and 
are useful in the discussion of persistence with parallel 
dynamics in other models. 
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